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Abstract 

The vacuum polarization of a quark, when considered in terms of the external mo- 
mentum q 2 , is a function of the Stieltjes type. Consequently, the mathematical theory of 
Pade Approximants assures that the full function, at any finite value of q 2 away from the 
physical cut, can be reconstructed from its low-energy power expansion around q 2 = 0. 
We illustrate this point by applying this theory to the vacuum polarization of a heavy 
quark and obtain the value of the constant K ^ governing the threshold expansion at 
order 0(a 2 s ). 



1 Introduction 



One interesting object to study in connection with the physics of heavy quarks is the vacuum 
polarization function of two electromagnetic currents. This requires high order perturbative 
calculations which, because of the obvious need to keep a nonzero mass m for the quark, 
become extremely difficult to perform. This is why, while the 0(a®) and 0{a\) contributions 
have been known for a long time pQ , state of the art calculations can only produce a result 
at 0(a 2 ) in the form of an expansion at low energies (i.e. q 2 = 0), at high energies (i.e. 
q 2 = oo) or at threshold (i.e. q 2 = 4m 2 ), but not for the complete function, which is still 
out of calculational reach. In this circumstances, it would of course be very interesting to 
be able to reconstruct this function by some kind of interpolation between the former three 
expansions. 

After the work in refs. 0E1I3], it has become customary to attempt this reconstruction of 
the vacuum polarization function with the help of Pade Approximants. Since these approxi- 
mants are ratios of two polynomials in the variable q 2 , they are very suitable for the matching 
onto the low-energy expansion. This is so because this low-q 2 expansion is truly an expansion 
in powers of q 2 , as a consequence of the finite energy threshold.] starting at 4m 2 . However, it 
is clear that they cannot fully recover the nonanalytic terms which appear, e.g., in the form 
of logarithms of q 2 in the expansion at high energies (or at threshold, where there is also a 
squared root behavior). Therefore, what is really done in practice is to first subtract all these 
logarithmic pieces from the full function (impossible to match exactly with a Pade) with the 
help of a guess function with the appropriate threshold and high-energy behavior, and then 
apply Pades to the remaining regular expressiorH. In this way, the authors of Ref. [5j were 
able to compute, e.g., the value of the constant appearing in the C(a 2 ) expansion of the 
vacuum polarization at threshold, which has not yet been possible to obtain from a Feynman 
diagram calculation. 

Although this result is very interesting, the construction is not unique. As recognized 
in Ref. [5], some amount of educated guesswork is required in order to resolve the inherent 
ambiguity in the procedure. For instance, a certain number of unphysical poles are encoun- 
tered, and some additional criteria have to be imposed in order to decide how to discard these 
poles. Since the resulting ambiguity leads to a systematic error which needs to be quantified, 
this error is then estimated by varying among several of the possible arbitrary choices in the 
construction. Although all these choices are made judiciously and in a physically motivated 
manner, it is very difficult to be confident of the error made in the result, which obviously 
has an impact on the value extracted for the constant K^ 2 \ 

In this short note we would like to point out that, regarding the vacuum polarization 
function, one can do away with all the above ambiguities. The vacuum polarization function 
belongs to a class of functions (the so-called Stieltjes functions) for which a well-known 
theorem assures that the diagonal (and paradiagonal) Pade Approximants converge. The 
result of this theorem together with the fantastic amount of information obtained on the 
Taylor expansion around q 2 = 0, for which 30 terms are known [6], will allow us to predict a 
value for K^ 2 ' . As it turns out, our result is very close to that of Ref. [5], although slightly 
smaller. 

The Theory of Pade Approximants is sufficiently developed to make these approximants 

1 In perturbation theory, this is true so long as purely gluonic intermediate states are not considered. Beyond 
perturbation theory, the threshold occurs at 4m£, where m n is the pion mass. 

2 In fact, this is done after a conformal mapping whereby all the (cut) complex plane is mapped into a circle 
of unit radius. 
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a systematic mathematical tool. As we will see in this article, and has been exploited in 
Refs. [21 El IU [5], they can be very useful for higher order calculations in perturbation theory. 
Furthermore, they can also be a very interesting conceptual tool for large- iV c QCD [7j. In 
this case this is true even when the function is not Stieltjes because large- N c QCD Green's 
functions are meromorphic and there are powerful theorems in this case as well [8]. Finally, 
Pades can also be used for analyzing the experimental data by fitting to a rational function 
rather than the more commonly used polynomial fitting [9|. They are also instrumental in 
discussing methods of unitarization |10j . 



2 Pades and the vacuum polarization 

Let us start by defining the vacuum polarization function H(q 2 ) through the correlator of two 
electromagnetic currents j^(x) = q(x)j fl q(x), 

{g^q 2 - q^u) I% 2 ) = -if dx e^ x ( \T j^(x)j u (0)\ ) , (1) 

where q^ is the external four-momentum. As it is known, due to the optical theorem, the 
e + e~ cross section is proportional to the imaginary part of II. As a result, Imll is a positive 
definite function, i.e. 

Imn(t + ze)>0, (2) 

a property which will become crucial in what follows. 

In perturbation theory II(g 2 ) may be decomposed to 0(a 2 ) as 

u( q 2 ) =n(% 2 ) + (^) n( 1 )(g 2 ) + (^) 2 n( 2 )(, 2 ) + otf). (3) 

For definiteness, a s denotes the strong coupling constant in the MS scheme at the scale 
/i = mp i e , but this is not important for the discussion which follows. Equation ([3]) will be 
understood in the on-shell normalization scheme where a subtraction at zero momentum has 
been made in such a way as to guarantee that 11(0) = 0. 

As it is well known, the vacuum polarization in Eq. ([3]) satisfies a once subtracted 
dispersion relation, i.e. 

2 f°° dt 1. 



n <<rt = '7 t(t _,._ ig) -i-m+fc). (4) 

Since all diagrams with intermediate gluon states are absent up to 0(a 2 ), the lower limit 
for the dispersive integral starts, in fact, at a finite value given by the threshold for 
pair production, i.e. 4m 2 . This fact only carries over to higher orders in a s provided these 
intermediate gluon states are neglected. From now on, we will restrict ourselves to the vacuum 
polarization in Eq. ([3]) to C(a 2 ), neglecting higher orders in a s . 
In terms of the more convenient variable 

z 55 h ' (5) 

one can rewrite Eq. (HJ), after redefining u = Am? /t, a^l 

p X j 

U(z) = z — : Imn(4m\ _1 + ie) . (6) 

Jq 1 — UZ — is TT 
3 We are simplifying the notation by replacing H(4m 2 z) — > 11(2;). 
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Recalling that a Stieltjes function is defined as jllpl 



/(*) 



l/R 



d(j)(u) 



uz 



(7) 



where 



is any nondecreasing function, one sees that the identification 



1 



u) = — Imll (4m 2 w 1 + ie) du 



(8) 



allows one to recognize that the integral in Eq. © defines the Stieltjes function z _1 n(z). 
As one can see, the positivity property Eq. ([2]) is crucial for the identification (|8|) to be 
possible. The representation of the function f(z) in Eq. (J7]) clearly shows a cut in the z 
complex plane on the positive real axis for R < z < oo. For the physical function II(z), this 
of course corresponds to the physical cut in momentum for 4m 2 < q 2 < oo, i.e. the physical 
case corresponds to R = 1 in Eqs. (|7|8|) . Furthermore, just like the function f{z) in Eq. ([7]) 
has a power series expansion convergent in the disk \z\ < R, so does the function IT(z) in Eq. 
^ have a power series expansion convergent in the disk \z\ < 1. 

A Pade Approximant to a function f(z), which will be denoted by Pff (z), is the ratio of 
two polynomials of degree M and N (respectively such that its expansion in powers of z 
about the origin matches the expansion of the original function up to and including the term 
of 0(z M+N ). When the original function f(z) is Stieltjes with a finite radius of convergence 
about the origin, R, it is a well-known result in the theory of Pade Approximants that the 
sequence P^j +J (z) (with J > —1) converges to the original function, as N —* oo, on any 
compact set in the complex plane, excluding the cut at R < z < oo [TT]. This excludes , in 
the physical case, the cut at 4m 2 < q 2 < oo (recall that R = 1). The position of the poles in 
the Pade Approximant accumulate on the positive real axis starting at threshold, q 2 = 4m 2 , 
mimicking the presence of the physical cut in the original function. When Pades are applied 
to the vacuum polarization, this means, in particular, that there can be no spurious pole 
outside of the positive real axis in the z plane and, consequently, no room for ambiguities. 
Furthermore, the convergence of the approximation (and the error) can be checked as a 
function of N, as we will see. 

3 Analysis 

In Eq. ©, the full functions Ii^ 0,1 ) (q 2 ) are known. They are given by the following expressions 



n<°>(z) 
nW( 2 ) 



16vr 2 
3 

167T 2 



20 4_ 4(1 -,)(! + 2z) 
9 3z 3z V 1 



5 + 13 _ (l-s)(3 + 2z) + (l-s)(l-16s) c2 

6 6z z " ' 6z 1 ' 



(l + 2z) ( , ,d \ I(z) 

t 1+23(1 -')s)-r 



(9) 



4 In Ref. [TT], the variable is chosen to be —z rather than z. 

5 Without loss of generality, the denominator polynomial of degree N is chosen to be unity at z — 0. 
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where 



Hz) 



G(z) 



6 



C 3 + 4Li 3 (-u) + 2Li 3 (w) -8 2Li 2 (-«) +Li 2 (u) 



lnu 



2 ln(l + u) + ln(l - it) 



In 2 tt , 



2u\au 
v?-l 



with 



u 



Vi-iA + i 



(10) 



However, the situation with the function n( 2 )(g 2 ) is different. In fact, Il( 2 )(g 2 ) is only 
partially known through its low-energy power series expansion around q 2 = 0, its high-energy 
expansion around q 2 = oo and its threshold expansion around q 2 = Am 2 , but the full function 
has not yet been computed. Unlike the latter two expansions, for which only a few terms are 
known, our knowledge of the expansion of H^ 2 \q 2 ) around q 2 = is very impressive, after 
the work of Ref. [6] where 30 terms of this expansion were computed. 

Although the full vacuum polarization function H(q 2 ) is Stieltjes, there is no reason why 
all the individual contributions n> 0,1 ' 2, "")(g 2 ) should also have this property. Amusingly, 
however, this happens to be true both for U^(q 2 ) and H^(q 2 ) pffi As we will now show, 
this is no longer the case for H^ 2 \q 2 ) because its power series expansion around q 2 = does 
not satisfy certain determinantal conditions which hold for a Stieltjes function. 

Defining the power expansion around z = of the Stieltjes function f(z) in Eq. ([7]) as 



f(z) = J2fnZ n 



(11) 



n=0 



the coefficients f n satisfy the following determinantal conditions. Let D(m, n) be the deter- 
minant constructed with the Taylor coefficients f n 



D(m, n) 



fm frn+1 
fm+1 fm+2 

fm+n /m+ra+l 



fm+n 
fm+2 n 



A Stieltjes function must satisfy D(m,n) > 0, for all m,n [11] ■ However, using the f n 
coefficients given in Ref. [H] (in the on-shell scheme, with the number of light flavors ri£ = 3): 



0.631107 + 0.616294 z + 0.56596 z 2 + 0.520623 z 3 + . . . 



(13) 



one can immediately see that, e.g. 



D(0 J 1) = L>(0,1) 



0.631107 
0.616294 



0.616294 
0.56596 



-0.0226376 < 



(14) 



This proves that the individual function 11^ (q 2 ) is, all by itself, not a Stieltjes function, 
even though the combination Ii(q 2 ) in ([3]) is. Therefore, we will now focus on applying the 
Theory of Pade Approximants to the full combination H(q 2 ) in Eq. ([3]) in order to extract 
information on the individual term H^ 2 \q 2 ). 



6 The case of IT^(q 2 ) is trivial as it coincides with the full vacuum polarization for 
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Figure 1: Number of decimal places reproduced by the rational approximation in Eq. H5\) to the 
function U.^'(z) as a function of z, in the interval 0.4 < z < 0.9. 



Since, as it is obvious from Eq. ([3]), the function H(q 2 ) depends on the value of a s , any 
Pade Approximant to it will also depend on the value of a s , i.e. Pj^ +J (z;a s ). This means 
that it is possible to construct a rational approximation to the three functions H^ ' 1 ' 2 )^ 2 ) 
from three different sequences of Pade Approximants to H(q 2 ) constructed at three arbitrary 
values of a s , let us say a s = 0, ±/3, with (5 sufficiently small so as to be able to neglect the 
terms of 0(af) in Eq. ([3]) . In this way one obtains 



n(°)(z) 
nW(z) 



p 



% +J (z;a s = 0) 



IT 



23 [ P^ J (z;a s =f3)-P^ +J {z-a s 

'2 



n(2)(z) w ^{p$ +J te<*. = fl+P% +J te<*. 



-0)} 



r,pN + J 



(z;a s = 0)} , (15) 



where J > — 1 and N — > oo. Since the value of (3 chosen is arbitrary, the N — > oo limit 
should produce results which are independent of /3, due to the convergence of the Pade 
Approximants to Yl(q 2 ). Therefore, one should see that the three combinations (|15p are 
increasingly independent of (5 as N growsJll This is indeed what happens. 

Furthermore, since we know the exact function IT 1 )^), we can compare it to the rational 
approximation on the right hand side of the second Eq. (I15D in order to test the approxima- 
tion. Figure [U shows the number of decimal places reproduced by this rational approximation 
in the interval 0.4 < z < 0.9, when = 14 and J = (i.e. the diagonal Pade P^), for 
values of (3 in the interval 0.1 < (3 < 1. As one can see, the dependence on [3 cannot be 
distinguished, and the accuracy reaches, e.g., ~ 10 decimal places at z = 0.9. 

The third Eq. (|15j) yields the desired approximation to H^ 2 \z). To be precise, it gives us 
a rational approximation to Ii^ 2 \z) in any compact set of the z complex plane, away from 
the cut 1 < z < oo. Since the threshold expansion at z ~ 1 can be written as p2 [12] 



1.72257 



+ [0.34375 - 0.0208333 n e ] ln 2 (l 



+ [0.0116822 n e + 1.64058] ln(l - z) + AT (2) 
+ [-0.721213 - 0.0972614 n e + 3.05433 ln(l - z)\ y/l-z + 0(1 - z) , (16) 

in terms of an unknown constant K^ 2 \ our Pade Approximation (|15j) may be used to deter- 
mine this constant, as we will next discuss. In this threshold expansion we take ri£ = 3 as 



7 This independence of @ in the case of H^^q 2 ) is trivially true. 
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Figure 2: Predicted value for K^ 1 from the sequence of diagonal, Pff (left panel), and first paradiag 
onal, f^ -1 (right panel), Fade Approximants. This figure corresponds to (3 = 0.5. 



the number of light flavors. Even though the numerical coefficients have been rounded off for 
simplicity, they may be extracted exactly from the results in Ref. |12j . 

Since Pades are not convergent on the physical cut, it is impossible to match the rational 
approximants (|15p to the threshold expansion (|16p as a function of z. This fact is obvious from 
the presence of logarithms and squared roots in Eq. (|16|) . However, both approximations 
(I15|) and (|16j) are valid for values of z at a finite distance from the cut and, in particular, in 
a certain window in the interval < z < 1. Within this window, a numerical matching of 
(|15l) and fjl6[> is possible and will in fact allow us to determine the unknown constant K^ 2 \ 

In order to determine this window, we make the following observations. First, although 
the rational approximation (|15p is convergent as N — > oo in the interval < z < I, it is more 
accurate the closer one gets to z = in this interval, for a given value of N. On the other 
hand, the threshold expansion f) 161) is more accurate the closer one gets to the branching 
point at z = 1. From these two competing effects it is possible to determine an optimal 
window in z by minimizing a combined error function. We will call this error function £{z). 

The function £{z) has to take into account the error from the Pades as well as the 
error from the threshold expansion. To estimate the error from the Pades, we consider the 
difference between two consecutive elements in the sequence, i.e. \P^ +J — P$Zi +J \- As to 
the threshold expansion, we estimate its error as |1 — z\, since the expression (|16p is accurate 
up to terms of 0(1 — z). Therefore, in order to avoid possible accidental cancelations between 
the two errors, we define our combined error function as the following sum: 

S(z)= ^{P% +J (z;a s = P)+P% +J (z;a s = -(3) - 2P% +J (z; a s = 0)} - {iV^iV-l} 

+ (17) 

Minimizing £ (z) with respect to z in the interval < z < 1, for every given values of N and 
f3, we may determine a value of z at the minimum, namely z*. This z* is then the one used 
to determine the constant as 



R{2) ^ W { P N +J( > Z *> as =® + P N +J ( z *-> <*s = ~P) ~ 2P% +J (z*; a s = 0)} - u[ 

(18) 

for the given N and (3. In Eq. (US), n^(z*) stands for the expression in Eq. ([16[) without 
the constant K^ 1 and, of course, without the term 0(1 — z), evaluated at z = z* . The 



6 




0.5 0.6 0.7 0.8 0.9 1.0 



Figure 3: Matching of the rational approximant in Eq. H5\) (solid blue line) to the threshold expansion 
^th * n (solid-dashed red line) for N = 14, J = and the value of in H9\). This figure 
shows the result for different values of (3 in the range < (3 < 1, but the dependence on (3 is so small 
that cannot be discerned. 

knowledge of 30 terms from the low-energy expansion gives us enough information to be able 
to construct up to the Pade Pjj* from the diagonal sequence, and up to the Pade P-J5 from 
the first paradiagonal sequence. This corresponds to J = and J = — 1 in Eq. (|18|) . In 
all cases considered we have varied (3 in the generous range < f3 < 1, but our results are 
insensitive to this variation within errors, as expected. 

4 Result 

Our results for the constant K^' are shown in Fig. [2. This figure shows the convergence of 
the diagonal sequence and the first paradiagonal sequence as a function of the order in the 
Pade. As one can see, we find a very nice convergence in both cases, with compatible results. 
Based on this analysis, we obtain the following value of A"' 2 ': 

= 3.71 ± 0.03 . (19) 

This result is very close to, although slightly smaller than, the value obtained in Ref. [5], i.e. 
= 3.81 ± 0.02. This is the main result of this work. 

The error bars shown on Fig. [5] have been calculated as ±£(z*). Looking at the figure 
we see that the change in the value of from one element of the sequence to the next is 
included in the errors shown, which we interpret as a sign that the estimate for the error we 
have made is rather accurate. Although we have taken symmetric errors for simplicity, it is 
also clear from the figure that the approach to the true value is made from below, so that a 
slightly more accurate determination could be achieved with the use of an asymmetric error. 
Apart from that, given the present knowledge of the expansions at low energy (|13p and at 
threshold (|16p . we find it difficult to believe any error estimate which could significantly go 
below our figure in Eq. (|19p . Of course, should more terms in either expansion be known, a 
rerun of our analysis could immediately produce a more precise determination of . 

Figure [3] shows the matching of the rational approximant (i.e. the right hand side of the 

_____ (2*\ _____ 

third of the Eqs. (|15p ) to the threshold expansion given by ITv (z) in Eq. (|16|) . for N = 14 
and J = 0, i.e. with the Pade P/4 , and for the value of we have obtained. As one can 
see, this Pade is able to reproduce, with high accuracy, the threshold expansion behavior in 
a window 0.92 < z < 1. At z = 1 and above, the two lines in Fig. __will again diverge from 
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Figure 4: Location of the poles in the pade in the complex plane. Notice the accumulation of poles 
at z > 1, simulating the physical cut. 



each other, just as they do at low z. The value of z* minimizing the error function £{z) in 
(I17|) was found at z* ~ 0.98 in this particular case. 

For illustration, in Fig. H]we show the position of the poles in the Pade P^. As one can 
see, all the poles are sitting on the positive real axis above z = 1, as it should be. Notice 
how they accumulate in the region z > 1. This is how Pades approximate the physical cut 
present in the original function. As ensured from Pade Theory, this behavior was found in 
all the Pades considered. 

Finally, as a further test of our method, we have calculated the value of the constants 
and h[ 2 ^ which appear in the large- z expansion of the function U^ 2 \z) (we take nj = 
ri£ + 1 = 4 in the following expression): 



U m g h~z( z ) = (0-034829 - 0.0021109 n f ) ln z (-4z) + (-0.050299 + 0.0029205 n f ) ln(-4z) 
+ H {2) + (0.18048 - 0.0063326 n/)^— ^ + (-0.59843 + 0.027441 n/)^— ^ 

tt{2) 

+ — ^ + 0^- 3 ln 3 (-z)J . (20) 

(2) 

Using our method, we find Hq ' = -0.582 ± 0.008. This result is to be compared to the 

(2) (2) 

true value Hq ' = -0.5857 @]. If we now input this exact value of .Hq , by a rerun of the 

(2) (2) 

method, we may then determine the value of . In this way, we findFf' = -0.194±0.033, 

(2) __ 

which is to be compared to the exact value H[ = —0.1872 [4j. Again, we find this agreement 
rather reassuring. 
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